Draft version February 4, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



o 
o 

(N 

X) 
D 



6 



> 

00 

o 

O 

o 



X 
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ABSTRACT 

We present the first high-resolution X-ray study of emission Hne variabihty with superorbital phase 
in the neutron star binary LMC X-4. Our analysis provides new evidence from X-ray spectroscopy 
confirming accretion disk precession as the origin of the superorbital period. The spectra, obtained 
with the Chandra High-Energy Transmission Grating Spectrometer (HETGS) and the XMM-Newton 
Reflection Grating Spectrometer (RGS), contain a number of emission features, including lines from 
hydrogen-like and helium-like species of N, O, Ne, and Fe, a narrow Ovii RRC, and fluorescent 
emission from cold Fe. We use the narrow RRC and the Hea triplets to constrain the temperature 
and density of the (photoionized) gas. By comparing spectra from different superorbital phases, we 
attempt to isolate the contributions to line emission from the accretion disk and the stellar wind. There 
is also evidence for highly ionized iron redshifted and blueshifted by ~ 25, 000 km s~^. We argue that 
this emission originates in the inner accretion disk, and show that the emission line properties in LMC 
X-4 are natural consequences of accretion disk precession. 

Subject headings: accretion, accretion discs — pulsars: individual (LMC X-4) — binaries (including 
multiple): close — stars: neutron — X-rays: binaries 



1. INTRODUCTION 

The eclipsing binary LMC X-4 consists of a 1.25 M© 
accretion-powered X-ray pulsar orbiting a 14.5 Mq 08 
HI star with a period of 1.4 days (Kelley et al. 1983; 
van der Meer et al. 2007). The system exhibits a 13.5- 
s pulsation due to the spin of the neutron star and a 
long-period superorbital variation of approximately 30 d 
(Lang et al. 1981; Ilovaisky et al. 1984), over which the 
X-ray flux may vary by as much as a factor of 60. Sim- 
ilar behavior has been observed in other X-ray binaries 
(Her X-1, Giacconi et al. 1973; SMC X-1, Gruber & 
Rothschild 1984). 

The standard model for the superorbital period in 
LMC X-4 is a precessing accr etion disk that periodically 
obscu res the neutron star (jHeemskerk fc van Paradija 
Il989l ). For high- luminosity systems, radiation pressure 
can generate a stable warp in the accret ion disk, whic h 
then precesses due to radiation torque l|Pringld ll996D . 
During the low state, an observer views the edge of the 
disk and the neutron star is obscured; during the high 
state, t he disk face is r evealed and the neutron star is 
visible (|Pettersonlll977f l. Signatures of this obscuration 
are well documented in timing (Her X-1, Petterson et 
al. 1991; SMC X-1, Hickox et al. 2005; Neilsen et 
al. 2004) and broadband spectral (SMC X-1, Vrtilek 
et al. 2001; LMC X-4, Woo et al. 1995) studies of pul- 
sar binaries. Furthermore, high-resolution X-ray studies 
of the accretion disk corona in Her X-1 indicate that 
the disk is in fact o bserved edge-on during the low state 
(|Jimenez-Garate"et al. 2005 ) • 
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There is also substantial evidence to suggest that emis- 
sion lines in Her X-1 and LMC X-4 vary with superorbital 
phase ^ (Jimenez-Garate et al. 2002, hereafter JG02; 
Naik & Paul 2003; Naik & Paul 2004). However, the 
relationship between lines and the continuum is com- 
plex and may require contributions from the accretion 
column, the accretion disk, a radiatively-driven stellar 
wind, or an X-ray heated wind from the accretion disk 
atmosphere. The extent to which these emission lines 
vary with 5" depends on the true physical origin of the 
superorbital period. 

High-resolution spectroscopy is necessary to decom- 
pose the spectrum of LMC X-4 into its physical com- 
ponents. Particularly useful is the Hea emission line 
triplet (from helium- like ions). This triplet consists of 
the closely-spaced resonance line (r: Is^ ^Sq - ls2p ^Pi), 
intercombination lines (i: Is^ ^Sq - ls2p '^P2.i), and 
the forbidden line (/: Is^ ^So - ls2s =^Si). The ratio 
R = f /i can be used to determine the electron density of 
the emitting plasma for 10^ cm~'^ < Ue < 10^^ cm~^, and 
G = {f + i)/r is sensitive to the electron temperature in 
the range 10*^ K < Te < 10^ K (Gabriel & Jordan 1969; 
Porquet & Dubau 2000). Furthermore, because photo- 
and coUisional excitation populate atomic levels differ- 
ently, these line ratios can be used to infer the dominant 
excitation mechanism in astrophysical plasmas. 

We have carried out detailed spectral analysis of the 
first high-resolution X-ray grating spe ctra of LMC X- 
4, ob tained with the Chandra HETGS (ICanizares et al 
20051) and the XMM-Newton RGS (den Herder et al 



20011 ). These observations probe the transition from the 



low state to the high state of the superorbital period, 
revealing the physical processes at work. We present 
the first study of the plasma conditions in LMC X-4 
with density and temperature diagnostics from helium- 
like species; we exploit the variability of these lines to 
constrain the geometry and dynamics of the line-emitting 
regions. In §2, we describe the observations and data 
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Superorbital Phase 

Fig. 1.— RXTE/ASM lightcurve of LMC X-4, folded on the 
long-term period of 30.3 d; two cycles are shown for clarity. 'I' = 1 
is defined to be the peak of the high state. The Chandra HETGS 
grating observations of the state transition are shown as horizontal 
arrows, and the XMM-Newton RGS observation of the high state 
is shown as a vertical arrow. 

reduction; in §3 we discuss continuum fitting and line 
searches. We present and interpret the emission hues 
from LMC X-4 in §4 and §5, and conclude in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

LMC X-4 was observed with the RGS on XMM-Newton 
on 2003 September 9 (10:11:34 UT) for 113.5 ks (ap- 
proximately one full orbital period). This observation 
took place during the high state of the 30 d superorbital 
cycle. The transition from low state to high state was 
observed with the Chandra HETGS on 2007 August 30 
(01:52:21 UT), August 31 (10:49:17 UT), and September 
1 (16:18:44 UT) for 47.09 ks, 50.96 ks, and 44.06 ks, re- 
spectively, for a total exposure time of 142.11 ks on the 
source. The observations are reported in Table 1 and are 
shown with respect to the superorbital period in Figure 
1. 

The RGS data were processed with SAS v. 7.1.0 and 
standard data reduction tools. We obtained spectra from 
RGS-1 and RGS-2 with FWHM spectral resolutions of 
roughly AA = 0.06 A. Bad columns, including CCD chip 
7 in RGS-1 and CCD chip 4 in RGS-2, were flagged and 
subsequently ignored. We also extracted the EPIC-PN 
spectrum in order to constrain the power law index above 
3 keV in the high state. The EPIC-PN spectrum wiU be 
analyzed in a future paper (K. Dennerl, in preparation). 
For reasons to be discussed in §4.1, we also extracted a 
separate RGS spectrum of LMC X-4 during eclipse. 

The HETGS data were analyzed with the CIAO anal- 
ysis suite, v. 3.4.3. After reprocessing and filtering the 
data, we extracted Medium Energy Grating (MEG) and 
High Energy Grating (HEG) spectra, which have FWHM 
spectral resolutions of 0.023 A and 0.012 A, respectively. 
Because the destreak tool can cause spectral artefacts for 
bright sources, we used the order-sorting routine to re- 
move the ACIS S-4 streak. To achieve sufficient S/N, 
we combined all three HETGS observations. This may 
involve averaging over real phenomena, but the third ob- 
servation is close to an order of magnitude brighter than 
the other two combined, so the effect should be small. 
In general, we do not account for orbital variability. We 
performed the following spectral analysis in ISIS (Houck 
& Denicola 2000; Houck 2002). 
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Fig. 2. — Data: model ratio for the RGS spectra near 12 A, 
showing the broad bump and the location of emission features from 
H-like and He-like Ne. The bump obscures the Neix lines and 
prevents a completely reliable determination of the Ne X line flux. 
In our continuum flts, we model it as a broad Gaussian with AA ^ 
1.3 A and equivalent width Wq ~ 60 eV. 

3. CONTINUUM SPECTROSCOPY AND LINE SEARCHES 

Historically, the X-ray continuum in LMC X-4 has 
been modeled as an absorbed hard power law plus a 
soft excess, which is evident because of the low Nh to- 
wards the LMC (|Hickox et al.ll2004D . The power law 
likely consists of Comptonized photons from the neutron 
star. Its slope varies between P ~ 0.5 — 0.7; the soft 
excess has been described as a combina tion of thermal 
Comptonization (ha Barber a et aLll2001h . blackbody ra- 
diation (.N aik & Paul 2002), and thermal bremsstrahlung 
(|Naik fc Pa ul 2004). We set the minimum Nn at the 
Gala ctic value, 5.78 x 10^° cm~^ (|Dickev fc LockmanI 
Il990l ) and fit the spectra with various combinations of 
these hard and soft components. 

However, we are unable to find a completely physical 
model for the high-state continuum due to a broad bump 
near 12 A (see Figure 2; also see JG02 for a discussion of 
possible origins for this feature in Her X-1). This bump 
is apparent in both RGS detectors. We model it as a 
Gaussian with width AA '--^ 1.3 A and equivalent width 
Wq 60 eV. For this reason we will refrain from drawing 
any significant physical inferences from our continuum 
model. In order to achieve xt ^ 3.4 (d.o.f. 565) in 
the RGS, it is necessary to include narrow Gaussians 
to model the spectral lines. We include ten Gaussian 
components: six for the helium-like triplets of nitrogen 
and oxygen, two for the Lya lines of nitrogen, oxygen, 
and two unidentified strong absorption lines. The bump 
prevents reliable analysis of the Ne IX triplet and the Ne X 
Lyman-a line. We let the fluxes of the Gaussian lines 
vary when fitting the continuum, then fit all the line 
parameters separately. 

Among those tested, the best spectral model for 
the high state data is thermal bremsstrahlung with 



0.455 



-1-0.022 
-0.001 



keV plus the broad Gaussian 



bump, the narrow lines, a blackbody with /cTbb — 
43.4lJ;^ eV and a power law with P = 0.813 ± 0.007 
(constrained by the EPIC spectrum above 3 keV), ab- 
sorbed by the Galactic neutral hydrogen column (A^h = 
5.78 X 10^" cm~^). The absorption lines, which may be 
interstellar in origin, are found at 23.54llg A and 
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TABLE 1 
Observations of LMC X-4 



Obs. 
ID 


X-ray 
State 


Grating 


Start 
Time'' 


T 

-■- exp 

(ks) 


Orbital 
Phase'' 


0142800101 
9573 
9574 
9571 


High 
Transition 
Transition 
Transition 


RGS 
HETGS 
HETGS 
HETGS 


52891.42 
54342.08 
54343.45 
54344.68 


113.5 
44.1 
51.0 
47.1 


0.25-1.19 
0.28-0.71 
0.25-0.62 
0.12-0.53 



^ MJD =JD-2,400,000.5^ Computed from the ephemcris of lLevine et al.l 
ll2OO0lV 



23.7l" y A, and have equivalent widths of -0.96 ± 1.4 eV 
and -0.2lto \l 

eV, respectively. This model results in 
xi = 1-55 (d.o.f. 537), which is not formally acceptable, 
but there is very little structure in the residuals and it is 
not our goal to model every spectral feature in the RGS. 

Because the Chandra and XMM-Newton observations 
probe different X-ray states of LMC X-4, we fit their 
continuum parameters separately. Unfortunately, there 
is insufhcient S/N in the MEG/HEG to constrain the 
column density independently, so in fitting the HETGS 
spectrum we fix A^h at the Galactic value (also measured 
by the RGS). To ensure that our characterization of the 
continuum is accurate given the low S/N, we group the 
data to a minimum of 40 counts per bin and include a 
Gaussian emission line to account for the strong Fe Ka 
hne at 6.4 keV. The best fit for the HETGS data from 
0.5-10 keV is bremsstrahlung with kT = 0.43 ± 0.03 keV 
plus a shallow power law with F — 0.40 ± 0.04. The iron 
Hne appears at 1.941 ± 0.002 A and has Wq = 130 ± 
30 eV; the model gives xt = 1-21 for 886 degrees of 
freedom. Assuming a distance of 50 kpc to LMC X-4, we 
measure unabsorbed continuum luminosities of (0.5-10 
keV) Lx = 1-2 X 10'^^ ergs s~^ in the state transition and 
ix = 1.0 X 10^® ergs s^^ in the high state. 

After finding a satisfactory continuum model for the 
HETGS, we proceed to model the spectral lines. We 
include Gaussian components for the the triplets from 
helium-like species of oxygen, neon, and iron as well 
as Lyman-a lines from hydrogen-like nitrogen, oxygen, 
neon, and iron (these were the strongest lines visible in 
our spectra), and find 90% confidence limits. Because 
of somewhat low S/N, we find it necessary to rebin the 
MEG spectrum to 0.092 A (4x worse than the nominal 
resolution) for A > 16 A; the rest of the data are binned 
to the nominal resolution. Our final fit gives xt — 1-08 
for 1758 degrees of freedom. We show the line parameters 
in Table 2. The RGS-1 and RGS-2 spectra are plotted 
in Figure 3, and the MEG spectrum is plotted in Figure 
4. The RGS eclipse spectrum is plotted in Figure 5. 

In order to confirm our line IDs, we perform simple 
checks with analytic photoionization models from xs- 
TAR, usin g a power law with F = 0.4 as our ionizing 
spectrum (jKallman fc McCravlll983 ). We compare our 
spectrum with the XSTAR predictions for other lines and 
adjust the model to fit for the appropriate line redshifts 
and intensities. In this way we are able to estimate the 
dynamics and ionization state of the emitting plasma. In 
the case of the narrow iron lines in the HETGS, we find 
evidence for lines separated from Fexxv by ^ 25,000 
km/s. In order to assess the possibility that these lines 
correspond to a physical Doppler shift, we search for sim- 
ilar Doppler shifts of our other lines as well. This process 
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Fig. 3.— RGS high-state spectrum of LMC X-4; RGS-1 is plotted 
in the top panel and RGS-2 is plotted in the center panel. Note the 
strong Ovill and Nvil Lyo Uncs at 18.97 A and 24.78 A, respec- 
tively. Gaps in the spectra correspond to ignored columns. The 
red, green, and blue models correspond to the continuum, Region 
I, and Region II spectra. We show the residuals for RGS-1 in the 
bottom panel. 

results in a total of 21 Gaussian emission lines (station- 
ary or Dopplcr-shifted) and two radiative recombination 
continua (RRC) in the HETGS. 

4. EMISSION LINES 

The line emission from LMC X-4 is dominated by 
hydrogen- and helium-like ions. In particular, the 
transition-phase HETGS spectrum and the high-state 
RGS spectrum contain narrow emission lines from Nvi, 
Nvii, Ovii, Oviii, Neix, and Nex. Because of the ab- 
sence of L-shell iron emission lines and the presence of 
RRCs from Nvi and Ovii (HETGS only), we conclude 
that these are recombination cascade emission lines. As 
can be seen from Table 2, these lines are all consistent 
with little or no Doppler velocity. In the helium-like 
triplets, the intercombination (i) line dominates. We 
will discuss line ratios in §4.2. In the transition phase 
we detect a strong fluorescent Fe Ka emission line and 
recombination lines from Fexxv and Fexxvi. 

A comparison of the transition-phase and high-state 
spectra yields an interesting discovery: emission lines 
from helium-like nitrogen, oxygen, and neon do not 
evolve in the same way as emission lines from the cor- 
responding hydrogen-like species. Our analysis leads us 
to the conclusion that at least three emission regions are 
necessary to explain the line properties in LMC X-4. In 
the next sections, we will summarize this argument and 
discuss each region individually. Because of chip defects 
in the RGS (20-24 A and 10.6-13.8 A for RGS-2 and 
RGS-1, respectively) and low S/N in the MEG, we wiU 
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TABLE 2 

X-iiA^" Emission Lixf.s in LAIC X- 1 





Ao 


A 










Line(s) 


(A) 


(A) 


(km s^-"-) 


Flux 


(eV) 


(km s~i) Grating 



Ne IX Hea 



Ovii RRC 
Ovil He/3 

Ovii Hea 



N VI RRC 
Nvi Hea 



13.417 (r) 
13.552 (i) 
13.698 (/) 

16.78 
18.63 

21.602 (r) 

21.802 (i) 

22.10 (/) 

22.46 
28.78 (r) 
29.08 (i) 
29.53 (/) 



13.560 



Region I: X-ray Irradiated Stellar Wind 

< 0.8 
1.7j 



+0.004 
0.001 



18.65 ± 0.04 
21 fi'i3+0°'"' 

21.805 ± 0.005 
21.841°:°* 



22.52 ± 0.08 
29.11 ± 0.02 



190] 

280^ 
440 



30 



580 
600 
+100 
870 



40 ± 70 

4~Q+500 



8io: 



230 ± 190 



7+0.6 
-0.5 

< 0.6 

< 1.8 

+ 2.1 
1.6 
+4.1 
2.8 

< 4.9 

13.31:1^ 

< 4.0 



+ 1.3 
-0.9 
+ 1.8 
1.5 

< 1.7 



1.0 
5.0 



3.4 
6.0 



1.4 
4.3 

5.3 



+3.3 

-3.2 



3.7 

< 7.6 

< 6.5 
19.4 



+3.8 
3.7 

< 6.3 



+ 1.7 
1.3 
+2.6 
2.1 
+3.6 
2.5 

< 0.4 

n 5+''-'^ 

9+°-3 

< 3.3 
0.3 ±0.2 

< 6.1 

< 0.2 
0.6 ±0.1 

< 0.2 



< 704 



700+"'"' 
'LI'J_34Q 

< 620 

< 400 

< 1030 



< 370 



Region II: Outer Accretion Disk 



HETG 
HETG 
HETG 

HETG 
HETG 

HETG 
RGS 

HETG 
RGS 

HETG 
RGS 

HETGS 
RGS 
RGS 
RGS 



Fexxvi Lya 


1.781 


1 7QK+0.008 
J-. IOJ_0.012 


7-1 r,+1280 
' ^^'-'-2100 


< 1.4 


< 33.2 




HETG 


Fe XXV Hea'' 


1.850 (r) 
1.857 (i) 
1.868 (/) 


1 869+°°°^ 


3050+1510 
1910+1™° 

-|on + 1480 


3.711;? 


86.5128.6 


3180ll?«° 


HETG 
HETG 
HETG 


Fe Ka 


1.938 


-^■''^-^-0.002 


470lf,° 


5.8ll:J 


1 '^c^ 0+25.3 
iOO.z_23.o 


1200l«° 


HETG 


Nex Lya 


12.135 


19 ir;o+0.002 


5801^° 


0.9 ±0.3 


3.311:? 


28011?° 


HETG 


Oviii Lya 


18.97 


18.97 ±0.05 


90+760 


Q q + 2.3 

•^••^-1.8 


^•^-2.2 


820l^g° 


HETG 






18.95 ± 0.02 


-250 ± 280 


45.6 ±3.7 


4.7 ±0.4 


227(1+320 


RGS 


Nvii Lya 


24.78 


24 70+0-05 


c-0+600 






< 330 


HETG 






24.79 ± 0.02 


110 ± 180 




3.1 ±0.3 


13701?°° 


RGS 






Region ill: Inner Accretion Disk 








Fexxv Hea'* 


1.850 (r) 
1.857 (i) 
1.868 (/) 


-, (,+0.006 
J.. I -'-U_Q 008 


-23560t?l°o 
-24690t?°°o 
-26410t??°0 


1-0-1.0 


-, +29.4 
4d.i_24 g 


< 4320 


HETG 
HETG 
HETG 


Fexxv Hea'* 


1.850 (r) 
1.857 (i) 
1.868 (/) 


9 no5+° °°^ 


23950^1^20 
22820tlii° 
21100+11°° 


1 .5+0.8 


00 Q+17.8 
•i3.»_15.6 


iqiO+2240 


HETG 
HETG 
HETG 


O vii He/3 


18.63 


20 61+° °'^ 


30220lllf 


0(5+3.1 


pr + 3.1 
C>.0__2.1 


< 930 


HETG 


Oviii Lya 


18.97 


17.36l°;°i 


-265701^^° 


< 1.8 


< 2.8 


< 5200 


HETG 






20.92t°;°i 


29190l?«°o 


Q 1+3.2 
•J. ^-2.2 


9 Q+3.0 


< 1100 


HETG 


Ovii Hea 


21.602 (r) 


iq Q1+0.06 
iy.Ml_0 Q3 


-24450l|« 


3.0 


< 3.3 




HETG 




21.802 (i) 


20 ,0+0.09 

ZU.1J_(, g2 


-23970ti5^° 




9 1+2.8 


< 3450 


HETG 


Nvii Lyo 


24.78 


22 qi+"-»2 


-2347011^° 




6 6+'5-'' 
o.o_4.n 


< 2000 


HETG 



Note. — Errors quoted arc 90% confidence ranges for a single parameter. Aq : rest wavelength; A : measured 
wavelength; Avghift ■ measured Doppler velocity; Flux: measured line flux in units of 10"'^ photons s~^ cm~^; Wo : 
line equivalent width; (Tv '■ measured line width. The high state was probed by the RGS, and the transition phase 
was observed with the HETGS (see Figure 1). Notice that the Hnes from Regions I and II are all consistent with a 
very small Doppler velocity (as long as Fexxv in the high state is dominated by the forbidden line — see §4.3). 
^ Because the Fexxv triplet is unresolved, we report the Fexxv line with three velocities, which correspond to the 
iron velocity assuming the line is dominated by the resonance line, the inter combination line, or the forbidden line, 
respectively. 
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Wavelength (A) 

Fig. 4. — MEG transition-phase equivalent of Figure 3 showing 
data, fit, and residuals; again the red, green, blue, and magenta 
models correspond to the continuum. Region I, Region II, and Re- 
gion III, respectively. We have labeled the strongest lines. 

focus on the O vii triplet when discussing hehum-hke 
species. Note also that the emission lines are grouped 
according to region in Table 2. 

4.1. Line evolution 

The most striking feature of the variable emission lines 
in LMC X-4 is that while the emission from Nvii and 
O VIII is highly correlated with the continuum flux, the 
emission from Nvi and Ovii is not. In particular, 
while the equivalent widths of the Lyman-a lines from 
hydrogen-like N and O appear to be independent of su- 
perorbital phase (at the 90% confidence level), the equiv- 
alent widths of lines from the corresponding helium-like 
species decrease substantially when the continuum flux 
increases. Interestingly, the flux in the helium-like lines 
appears to be constant over superorbital phase. 

For example, O viii Lya has Wq 
sured by the HETGS in the transition and Wq = 4.7±0.4 
eV as detected by the RGS in the high state. In compar- 
ison, the equivalent width of the O vii i line decreases by 
a factor of 13, from Wq ~ 11.5^3 6 in the transition 
(HETGS) to Wo = O.dtoi eV in the high state (RGS). 
The flux in this line, however, does not change: we mea- 



4.1^2:2 eV as mea- 



sure (in units of 10 photons s cm ) / 



13.3^^2 

in the transition (HETGS) and / = U.GtH the high 
state (RGS). 

Further comparison of the transition state spectra to 
the high state spectra reveals a noticeable increase in 
the broadening of these hydrogen-like Lya lines relative 
to their helium-like counterparts: the helium- like triplets 
have no significant velocity width in either the transition 
or the high state. In the high state, the Nvil and Oviii 
Lyman-a lines have = 1370 km and 2270 km s~^, 
respectively. Because the RGS is able to constrain the 
width of narrower lines (Nvii is narrower than 370 km 
in the high state), we believe that this broadening is 
real. If it is thermal in origin, the implied temperature 
is T ~ 10® K; Doppler broadening arising from bulk mo- 
tion is more likely. These values should be compared to 
the transition-phase widths, which are < 330 km s^^ for 
Nvii and 820 km s^^ for Oviii. 

These results are difficult to reconcile with a single 




■NlS 



Wavelength (A) 



Fig. 5. — RGS-1 eclipse spectrum of LMC X-4. During the 
eclipse, the Ovill Lyo line weakens relative to the Ovil triplet, 
indicating that Region II is not visible during eclipse. Given that 
Region II is highly ionized, exhibits somewhat broadened emission 
lines, and is obscured during eclipse, we argue that it is located in 
the accretion disk. 

emission region with a variable ionization state. The fact 
that the hydrogen-like lines brighten by the same fac- 
tor as the continuum while the helium-like lines do not 
change suggests that the hydrogen-like lines and helium- 
like lines are produced independently, probably in two 
distinct emission regions with different ionization states. 
The fact that the O viii line weakens relative to the O vii 
triplet during eclipse (compare Figures 3 and 5) also sup- 
ports at least two emission regions. The presence of high- 
velocity (- 25,000 km s^^) iron during the state tran- 
sition implies that a third region may be necessary to 
explain all the lines in LMG X-4. 

The first region (hereafter Region I), is defined as the 
source of the helium-like triplets and weak hydrogen-like 
lines (note that this refers only to N, O, and Ne). The 
luminosity of Region I is not affected by superorbital vari- 
ations in the X-ray flux; if it consists of a photoionized 
plasma, this might indicate that its luminosity is tied to 
the intrinsic luminosity of the neutron star. Region II 
is the source of the strong Lya lines, whose intensity is 
tied to the observed X-ray continuum. This region makes 
little or no contribution to the helium-like triplets. Re- 
gion HI, for reference, is the origin of the Doppler-shifted 
Fexxv fines. 

4.2. Region I: the X-ray irradiated stellar wind? 

As discussed in the preceding section, we define Region 

1 as the origin of the helium-like emission triplets. These 
triplets are useful for plasma diagnostics because their R 
and G line ratios can constrain the density, temperature, 
and dominant ionization process in the emitting plasma. 
Given the striking differences between line evolution in 
Regions I and II, these quantities are particularly impor- 
tant for a complete description of this high-mass binary. 

In Nvi, Ovii, and Neix, we observe prominent in- 
tercombination lines, weak-to-moderate resonance lines, 
and weak forbidden lines; for all these species, R ^ 
(with upper limits of ^ 0.5 at 90% confidence). Colli- 
sions can produce this line ratio by exciting the 2 ^Si 

2 transition, effectively converting forbidden line 
photons into intcrcombination line photons. Thus if coU- 
sional ionization dominates, it is possible to find lower 
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limits on n using the the upper limits on R from the 
values in Table 2 and the relations h ip be tween R and 
n determined by iPorguet fc Dubaul (|2000f ). With this 
method, i? ~ implies a high density plasma in Region 
I. 

However, photoexcitation by the intense UV field of 
the secondary c an compete with collision al excitation to 
produce i? ~ ()Mewe fc Schriiveij[l978l ). The photoex- 
citation rate is given by 

7re 

Wf-,i = F^f^Josc, (1) 

mc ' 

where m and e are the mass and charge of the electron, 
c is the speed of light, /osc is the oscillator strength, and 
is the flux (in units of photons s~^ cm~^ Hz~^) at 
the wavelength of the 2 ^Si 2 ^-Pi,2 transition (1637 
A for O Vli). Following JG02, we assume a UV point 
source (so that the radiation field and photoexcitation 
rate Wf^i decrease with the square of the distance d 
from the secondary) and calculate the distance dcrit at 
which the photoexcitation rate equals the radiative decay 
rate Wf. Inside dcrit, photoexcitation destroys forbidden- 
line photons, so the density limits from i? ^ are not 
applicable. 

At orbital phase 0.74 (similar to our observations) and 
superorbital phase between 0.8 and 1.1 (i.e. in the low- 
to-high state transition or in the high state), Preciado, 
Boroson, & Vrtilek (2002) measured the UV flux at 1637 
A to be Fxf^^ = 2.2 ± 0.3 x lO^^^ ^j-gg g-i ^1^-2 ^-1 
Taking the radiative decay rat e and osci llator strength 
for O VII from .Drakc (197l|) andlCann&Jiakkar (1992), 
respectively, we calculate dcrit ~ 3 x 10^^ cm. This is 
mu ch larger than the orbital separation a = 9.3 x 10^^ 
cm (|van der Meer et al]l2007l ) , so photoexcitation cannot 
be ignored in explaining i? ~ 0. If we neglect collisions 
altogether, so that i? ~ is produced by photoexcitation 
alone, we find that Ri < dcrit (we can ignore the orbital 
separation because it is much smaller than dcrit)- We 
conclude that R cannot be used to constrain the density 
of Region I. 

It should also be possible to constrain the density using 
the G ratio, because at high densities the resonance line 
is enhanced by collisional depopulation of the ls2s ^Sq 
level. This effect results in a measureable decrease in G 
for n > 10^2 cm-3 for Ovii (|Gabriel fc Jordan! fl97l . 
We attempt to account for collisions by modeling the 
spectrum with XSTAR. We fit the Ovii triplet with a 
grid of models in n and ^, covering 1 < log^ < 2 and 
log(n/cm-3) ^ {11^ 12, 13, 14, 15, 16}, but find that the 
relatively low S/N of the HETGS spectra does not permit 
a statistically robust constraint on n. 

G can also be used to determine the ionization mech- 
anism in Region I. In the low-density limit, G = 4 
is the cutoff between photoionized and collisional plas- 
mas; G > 4 i ndicates that photoionization dominates 
(jPorguet fc Du bau 2000) . This cutoff is somewhat lower 
at high density. Because the r line is a combination of re- 
combination emission, resonant scattering into the line of 
sight, and resonant absorption, G should be interpreted 
with care (jWoidowski et al.ll2003D . This is particularly 
true during the state transition (observed with Chan- 
dra), when the pulsar is at least partially obscured and 
resonant scattering into the line of sight is almost cer- 
tainly important. In any case, the observed resonance 



lines in LMC X-4 are weak, and we are generally only 
able to place a lower bound on G. In the high state, we 
measure G > 2.0 using N and O; in the transition, we 
find G = 2.3li j and G > 1.5 for O Vii and Neix, respec- 
tively. These values allow for photoionization dominance, 
but we are are unable to use G alone to determine the 
dominant ionization process in LMC X-4. 

However, the comparison of hydrodynamic disk and 
wind models against UV spectra has shown that pho- 
toionization of the stellar wind can explain the orbital 
variability of emission and absorption line profiles in 
LMC X-4 (Vrtilek et al. 1997; Boroson et al. 1999; 
Boroson et al. 2001). Furthermore, we detect narrow 
radiative recombination continua (RRC) from Nvi and 
Ovii. In the narrow limit, the width of the RRC is pro- 
portional to the electron temperature. A joint fit to these 
two features indicates /cT ~ 1 eV, or T 1.2 x 10** K. 
The RRC are too faint to provide confidence limits on T 
with the redge model, but Gaussian fits (which are ac- 
ceptable because the features are narrow) imply cr < 4.3 
eV, or T < 5 X 10* K. Because collisional ionization is 
negligible at these temperatures, photoionization must 
be the dominant process in Region I. For reference, the 
continuum emission from Region I falls below the spec- 
tral sensitivity of the HETGS and RGS. 

Given these results, we suggest that the X-ray irradi- 
ated stellar wind is a promising candidate for the identity 
of Region I. We have already mentioned that this com- 
ponent contributes significantly to the UV spectrum of 
LMC X-4, and we have seen that the temperature of Re- 
gion I is typical of an 0-type star. At this point, we 
cannot legitimately distinguish the stellar wind from the 
face of the secondary. However, the simplest explanation 
for the apparent absence of variability on the superor- 
bital period in Region I is that it occupies a large solid 
angle as seen by the neutron star, and the stellar wind 
meets this criterion easily. We shall return to this topic 
in §5. 

4.3. Region II: the outer accretion disk/stream? 

Before discussing the Lyman-a lines from N vii, O viii, 
and Nex, it will prove useful to describe the iron lines 
in LMC X-4. During the state transition, we detect a 
number of iron emission lines. They are shown in Figure 
6, with Gaussian parameters in Table 2. It is straightfor- 
ward to identify three of these lines (at 6.4 keV, 6.7 keV, 
and 6.9 keV) as the Fe Ka complex, Fexxv Hea, and 
Fexxvi Lya, respectively. If we assume that the Fexxv 
emission is dominated by the forbidden line, then these 
three iron lines have the same velocity at 90% confidence, 
i.e. V km s^^. This is a reasonable assumption be- 
cause the critical density for Fexxv is well above 10^^ 
cm~"^, so its forbid den line will be strong fo r a very wide 
range of densities ([Gabriel fc Jordanlll972D . Two other 
lines at 6.1 keV and 7.3 keV constitute Region HI and 
will be addressed shortly. 

XSTAR photoionization models of this spectral region 
which confirm our identification of Fexxv Hea and 
Fexxvi Lya (see §3) indicate log^ ^ 3. At this ion- 
ization state, the contribution from helium-like species 
of N, O, and Ne is negligible. For the state transition 
in LMC X-4, XSTAR predicts Lya fines from N and O 
that are weaker than the Poisson noise in our spectra. 
Thus the weak narrow Lya lines that we detect during 
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Wavelength (A) 

Fig. 6. — HEG transition-phase spectrum of the LMC X-4 iron 
region. The red line is the continuum fit; blue and magenta corre- 
spond to Regions II and III, respectively. 



the state transition probably originate in Region I; their 
strong, broad counterparts in the high state originate in 
Region II. Disentanghng the exact contributions to each 
Hne from these two regions is beyond the scope of this 
paper; full photoionization modeling will follow in future 
work. 

At this point, we cannot proceed without associating 
the transition-phase iron lines with the high-state Lya 
emission lines from N, O and Ne. This association ap- 
pears to be confirmed by the presence of the Fe Ka and 
Fexxv lines in the EPIC high state spectrum (K. Den- 
nerl, in preparation). Furthermore, the iron line model 
easily provides the high ionization parameter required by 
the absence of helium-like species in Region II. In addi- 
tion, the comparable widths of the iron lines and the Lya 
lines, along with their comparable small Doppler veloci- 
ties, are consistent with a common origin. Thus we shall 
treat all these lines as originating in Region II. Note that 
the Lya lines become visible in the high state because of 
the correlation between the continuum flux and the flux 
from Region II. 

Our association of hydrogen- like N, O, and Ne with a 
hotter region of the binary also explains the absence of 
radiative recombination continua from these species dur- 
ing the high state. At an ionization parameter log^ ~ 3, 
the temperature in Region II is T > 5 x 10^ K {kT > 
50 eV). The corresponding RRC would be very diffi- 
cult to detect over the bright continuum of the high 
state. If the RRC norm alizations follow the scaling of 
iLiedahl fc Paerelj ()1996f ). then the absence of a strong 
Ovill RRC implies T > 10 eV at 90% confidence. Dur- 
ing the state transition, there is some evidence for a pos- 
sible narrow O viii RRC with kT = 0.5 eV (T - 6000 
K). Given the above discussion, we suggest that this fea- 
ture, if real, actually reflects the contribution of Region 
I to the Lya lines. 

Finally, we would like to determine the distance of Re- 
gion II from the neutron star. We have an estimate for 
the ionization parameter in this region from the iron line 
model, but no density estimate. However, the fact that 
Region II is obscured during eclipse implies a distance 
from the neutron star of Ru < i?*, where Rj, = 7.8 Rp) 
is the radius of the companion star (|van der Meer et al.l 
I2OOI . This is a weak constraint, but if we assume that 



the line broadening in Region II is due to orbital motion 
with a Keplerian velocity profile, then i?n ^ 1-5 x 10^^ 
cm, i.e. well inside the Roche lobe of the neutron star. 
We therefore conclude that Region II lies in the accre- 
tion disk. It should be pointed out that if the disk is 
cool and optically thick in the center but has a hot, opti- 
cally thin surface layer, then it would be possible for both 
Fe XXV and cold iron (for the Fe Ka line) to be present in 
Region II. This would also explain their similar velocity 
dispersions and radial velocities but different ionization 
states. However, if the width of the Fe Ka line is due to 
a blend of several charge states, then the cold iron might 
lie outside of Region II. 

4.4. Region III: the inner accretion disk? 

In addition to the three iron lines from Region II, we 
also detected emission lines at 6.1 keV and 7.3 keV; the 
XSTAR model for Region II reveals no obvious candidate 
IDs at the velocity of the three strong lines. However, 
the average of the line energies coincides with the Fe xxv 
Hea triplet, so we identify these lines as symmetrically 
relativistically Doppler-split Fexxv emission (with ve- 
locities of about 25,000 km s~^). The existence of a num- 
ber of faint unexplained lines near other strong emission 
lines detected at low energy with comparable velocity 
shifts (see Table 2) strengthens this interpretation. 

Similar Doppler splitting was observed by EXOSAT 
during a bright flare: pulse-phased spectroscopy of the 
iron line suggested the existence of two line components 
whose separation (^ 0.2c) varied with pulse phase and in- 
creased with the luminosity of the source (Dennerl 1989). 
The natural implication of that result is that those lines 
originated near the neutron star, possibly in the rapidly 
rotating accretion column. This would explain why the 
Doppler shift varied with pulse phase, and it is consis- 
tent with the visibility of the neutron star in the high 
state. In any case, it seems likely that this motion 
was dominated by the p ulsar's ~ lO^'^ G magnetic field 
(|La Barbera et al.ll2001l ). 

We consider the possibility that the Doppler splitting 
in the state transition is also caused by the pulsar's mag- 
netic field. Because the iron lines in the transition phase 
have velocity widths of order 1,000 km s~^, we can con- 
clude that they do not originate in the accretion column 
(otherwise they would be broadened by ^ 0.2c by the 
spin of the pulsar). However, material corotating with 
the magnetic field, e.g. at the inner edge of the accretion 
disk, could certainly be responsible for the observed lines. 
It is straightforward to calculate that material corotating 
with the neutron star at a speed of 25,000 km s~^ would 
lie at r > 5 X 10^ cm. This result could place Region III 
well inside the accretion disk. 

5. DISCUSSION 

In this paper, we identify a number of recombination 
emission lines in LMC X-4; by comparing observations 
at two epochs, we determine the evolution of these lines 
with ^' and find evidence for three emitting regions; here 
we discuss the implications of these results for the dy- 
namics and geometry of this binary system. 

5.1. Line variability and emission regions 

Our analysis of the variation of emission lines reveals 
clear differences in the properties of helium-like and 
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hydrogen-like N, O, and Ne in LMC X-4, which we de- 
tailed in the previous sections. We find it useful to sup- 
pose that these lines come from two different physical 
regions; we will argue that the disk precession model for 
the superorbital period, which reproduces the long-term 
lightcurve, predicts the existence of these very regions 
and their variability. This argument is depicted schemat- 
ically in Figures 7 and 8, which show the low-state and 
high-state geometry of LMC X-4 as seen from the orbital 
plane. 

First, we attempt to place constraints on the physical 
location of Region I. Without a robust density estimate, 
we cannot use the ionization parameter to estimate 
Nevertheless, in §4.2 we claimed that we can reasonably 
identify Region I with the X-ray irradiated face of the 
companion star or its substantial wind based on its vari- 
ability; we expand on this argument here. In order to 
conclude that Region I must occupy a large solid angle as 
seen by the neutron star, we must show that shadowing 
by the precessing disk is unlikely to produce a constant 
ionization parameter and line flux over the superorbital 
phase. 

To do so, we suppose that shadowing is important, i.e. 
that the illumination of Region I varies like the super- 
orbital lightcurve. Between the HETGS and RGS ob- 
servations, the number of ionizing photons in Region I 
increases by a factor of ~ 10, but its ionization parame- 
ter ^ remains constant. The relation 



implies that ni?^ must increase by a factor of 10 as well. 
But given the increase in the number of ionizing photons, 
Region I must also be smaller during the high state in 
order to produce the same line flux (in the optically thin 
limit). That is, the neutron star would have to illuminate 
a more dense, distant, and smaller region of the binary 
during the high state than it does during the transition. 

Clearly shadowing results in a very unusual picture of 
Region I. The accretion geometry is vastly simplified if 
shadowing is unimportant and Region I is constantly ion- 
ized by the intrinsic luminosity of the n eutron star, i.e. 
2.3 X 10^8 ergs s'^ (|Levine et al.ll2000D . In this case. 
Region I is very likely the X-ray irradiated wind from 
the companion. This explanation is consistent with the 
low temperature of the emitting plasma as well as the 
small line velocities in Region I and the importance of 
photoionization in LMC X-4. Furthermore, a relatively 
steady wind could provide a solid angle which does not 
change over time, so that its ionization parameter and 
line luminosity in Region I are independent of ^E". At 
present, it is difhcult to distinguish between the surface 
and the wind of the companion. But the precession of 
the disk should cause periodic illumination of the face of 
the secondary, so a more complete high-resolution study 
of line emission with ^! might be able to resolve this de- 
generacy and possibly constrain the inclination of the 
disk. 

While Region I is independent of wc find that Re- 
gion II, where the hot iron lines and the Lya lines origi- 
nate, is tied to the observed X-ray continuum. We argued 
in §4.3 that i?n ^ 1-5 x 10^^ cm and that Region II prob- 
ably lies in the outer accretion disk, which then fills ap- 
proximately 50% of the Roche lobe of the neutron star. 



Here it is relatively easy to suggest how the observed 
equivalent width changes might occur. In the outer disk, 
the neutron star is always visible, so the line luminosity 
is constant in time (just like in Region I). However, the 
projected area of the disk changes periodically, leading 
to stronger lines when more of the disk is visible, i.e. the 
high state. It is also possible that the lines originate in 
the accretion stream, which is periodically illuminated 
by the neutron star so that its line luminosity is tied to 
the observed X-ray flux. This effect results in stronger 
lines during the high state of the superorbital period. 

JG02 found evidence for a similar increase in emission 
line strengths from the low state to the short-on state 
in Her X-1. They suggested that the variation in line 
strengths might be explained by the enhanced visibility 
of the disk atmosphere during the short-on. While it ap- 
pears that Region II is a common feature of systems with 
precessing accretion disks. Region I seems to be absent 
in Her X-1. In fact, we can use this point as further 
evidence that Region I is the illuminated stellar wind. 
Her X-1 and LMC X-4 both have very similar orbital pa- 
rameters; the main dynamical difference between these 
systems is the mass of the companion star (~ 2.2 Mq in 
Her X-1 and 14.5 Mq in LMC X-4). Thus we expect a 
stellar wind to contribute much more to the spectrum in 
LMC X-4; the absence of a Region I equivalent in Her 
X-1 is not surprising if our hypothesis is correct. 

It seems that for a system with strong winds, the exis- 
tence of two emitting regions is a natural consequence of 
a tilted, precessing accretion disk, which illuminates dif- 
ferent components of the binary in different ways. One 
region presents a roughly constant solid angle for illumi- 
nation by the neutron star, so the resulting emission lines 
do not vary with 4*. A second component might precess 
with the disk or see the neutron star obscured, just as 
we do. By assuming the existence of these regions, we 
can sensibly decompose our observations into a physical 
picture of an accreting neutron star with a precessing 
accretion disk. 

5.2. Doppler shifts and the superorbital phase 

We also present evidence for the relativistic Doppler- 
splitting of the Fe xxv Hea emission line during the state 
transition in LMC X-4, and speculate that the lines orig- 
inate in the accretion disk, where highly ionized plasma 
may be dragged along with the pulsar's rotating mag- 
netic field at velocities of ^ 25,000 km s~^. We find 
evidence of comparably Doppler-shifted emission lines at 
low energy which are present during the state transition 
but absent in the high state. We consider here the ori- 
gin of such Doppler shifts and the implications of the 
Doppler evolution for the geometry of the binary. 

The conclusion that the Doppler splitting occurs in 
the accretion disk is highly dependent on the velocity 
field of the disk. For a Keplerian disk, orbital veloci- 
ties reach v ^ 25,000 km s~^ at tk S 1-9 ^ 10^ cm. 
But the disk around a highly magnetized neutron star is 
truncated at the radius of the magnetosphere Rm, where 
magnetic pressure balances the ram pressure of the gas 
(jFrank et al.ll2002D : 

Rm = 1-5 X 10Vi/^i?6 ^/■^Lg/ZVao^ cm; (2) 
m is the neutron star mass in solar masses, Rq its radius 
in 10^ cm, L37 its luminosity in 10'^^ ergs s~^, and /iao 
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Fig. 7. — A cartoon of the X-ray low state of LMC X-4. The disk is viewed roughly edge-on, obscuring the neutron star. A warp 
protrudes from the inner disk; as the ionized surface of this warp is dragged along by the pulsar's rotating magnetic field, we observe 
Doppler splitting of the disk emission lines. The apparent area of the outer disk is small, so lines from Region II are weak. Note that 
shadowing is not particularly important in Region I. This figure should be compared to Figure 8, which shows our high-state scenario. 



its magnetic moment in lO'^'^ G cm'^. With m = 1.25, 
Re = 1, Lsj = 23, and /xso = 10 (from lLa Barbera et al.l 
1200 It ), we find Rm ^ 2.3 x 10^ cm. Since the inner edge 
of the disk is much larger than rx, Keplerian motion is 
unhkely to produce the observed Doppler shifts. 

By assuming that Region III is corotating with the 
pulsar, we calculate Rm ~ 5 x 10^ cm. This is much 
larger than Rm, so the lines cannot originate at the very 
inner edge of the disk. However, i?iii might mark the 
edge of the inner disk warp, and a highly ionized sur- 
face layer of the disk might continue to corotate with the 
magnetic field out to this radius. Thus we conclude that 
the Doppler-shifted Region III lies in the inner accretion 
disk around LMC X-4. 

To clarify this model for visualization purposes, we 
consider a toy accretion disk, also shown in Figures 7 
and 8. Note that in these figures, the system is viewed 
from the or bital plane, not at the orb ital inclination of 
~ 68° ± 4° ()van der Meer et al.ll2007t ). This disk has a 
small but nonzero inclination with respect to the orbital 
plane, and precesses around the orbital angular momen- 
tum vector of the system. The inner disk is warped and 
slightly tilted with respect to the disk, so that the edge 
of the warp protrudes from the disk. The surface layer of 



this warp is highly ionized and corotates with the mag- 
netic field. Furthermore, because only a small part of 
this warp is visible the Doppler-shifted lines are narrow, 
not broad, relative to the line splitting. 

But we have observed LMC X-4 at multiple epochs, 
and it seems that the Doppler-shifted lines from Region 
III may not be visible at all superorbital phases. In fact, 
this Doppler evolution is a direct consequence of the pre- 
cession of the accretion disk: the changing orientation 
of the disk implies a maximum radial velocity in the ac- 
cretion disk during the low state, and a minimum in the 
high state. Thus Doppler-split lines should appear and 
disappear at different phases. It seems that we have the 
first dynamical confirmation of accretion disk precession 
as the origin of the superorbital phenomenon. 

Admittedly, the non-detection of Doppler-shifted lines 
in the high state RGS observation does not prove that 
the disk precessed from edge-on to face-on. As such, we 
address the possibilities that the Doppler splitting is not 
real (perhaps an effect of low S/N), or that we observed 
two distinct phenomena (and not emission lines evolving 
with ^I*). The Doppler-shifted Fexxv lines are both de- 
tected at > 90% confidence, so the conclusion is robust; 
the concern over S/N only applies to emission from low-Z 
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Fig. 8. — A cartoon of the X-ray high state of LMC X-4 to be compared with Figure 7 (the disk is more face-on now, and the pulsar is 
visible). Because the inclination of the disk has changed, the Doppler shifts from the inner disk/Region III are more difficult to observe. 
In contrast, the projected area of the disk is much larger, so we observe brighter lines from Region II. The solid angle of the wind seen by 
the neutron star does not change, so the line luminosity of Region I is constant. The hole in the inner disk represents the magnetosphere, 
where the magnetic field completely dominates the accretion dynamics. 

Fe Ka line, the Fexxv line, and the Fexxvi line are vis- 
ible from Region II. As the disk turns its edge to us, the 
continuum flux decreases and the Doppler lines appear, 
boosting the equivalent width of the iron li ne region. 

Thi s effect was already observed by iNaik fc Paull 
(|2003[ ) with the lower spectral resolution of the RXTE 
PCA. They found that the equivalent width of the iron 
complex was GSOlgg eV in the low state, 320 ± 30 eV in 
the transition, and 280 ± 20 eV in the high state. Our 
spectra are consistent with this result: in the transition, 
we find the combined equivalent width of the Fc Ka line 
and the Fexxv line to be 220 ± 40 eV; including the 
Fexxvi line, the equivalent width is Wq ^ 240 eV. If 
we include the Doppler lines, Wq increases to 300 ± 50 
eV (~ 310 eV when Fexxvi is accounted for). Thus the 
RXTE data are consistent with our claim that Doppler- 
shifted lines are undetectable in the high state, appear in 
the state transition, and peak during the low state. High- 
resolution observations of the low state are necessary for 
a definitive test of this interpretation. 

It may be that we observed two independent events 
in the state transition and the high state, but the ob- 
servations can be explained with a single model. The 



species. It is tempting to suggest that because the fluxes 
of the Lyman-a lines of hydrogen- like N, O, and Ne are 
tied to the continuum, the Doppler-shifted lines should 
be most detectable in the high state as well, but there 
are two reasons why this is not the case. 

First, the projected velocity of the accretion disk is 
smaller in the high state; if the disk inclination decreases 
sufficiently, this effect could make the Doppler-shifted 
lines difficult to distinguish from the stronger v ^ line. 
By fitting the RXTE ASM lig h tcurv e and P Cygni pro- 
files in the UV, 'B oroson et al.l ()1999l ) determined a disk 
tilt of 31° ± 2° off the orbital plane. Combined with the 
orbital inclination of 68°, this result implies that the ap- 
parent inclination of the disk warp varies from 100° 
to ^ 35° over the superorbital phase, so that Doppler 
velocities in the high state should be a factor of 2 
smaller than in the transition. Since our grating spectra 
can easily resolve such velocities, the projected velocity 
of the disk alone cannot explain the disappearance of the 
Doppler lines in the high state. 

Second, and more importantly, the equivalent widths 
of the Doppler-shifted lines should actually be anticorre- 
lated with the continuum flux. In the high state, only the 
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spectroscopic interpretation of the disk precession model 
completely explains the variability of the emission lines 
in LMC X-4 with superorbital phase; all we require is 
that the equivalent widths of the Region II lines be in- 
dependent of 5". 

This simple requirement can be met if Region II is in 
the accretion disk or otherwise very close to the neu- 
tron star. Then the observed flux in these lines will be 
strongly correlated with the X-ray continuum, and their 
equivalent widths will vary slowly with 4". Then as the 
disk precesses out of the high state, we observe addi- 
tional Doppler-shifted emission lines whose velocity and 
flux increase into the low state before disappearing again. 

6. CONCLUSION 

Only a handful of XRBs are known to have persistent 
superorbital periods, the most famous of which are He r 
X-1, SS 433, SMC X-1, and LMC X-4 (|Wen et al.ll200l . 
To date, there has been very little systematic study of 
the superorbital variation of the strengths or positions of 
X-ray emission lines from the precessing accretion disk in 
any of these systems (grating s tudies of SS 433, fo r exam- 
ple, have focused on jet knots (jLopez et al.ll2006f l). Such 
studies can be extremely enlightening, given the power 
of high-resolution spectroscopy to reveal the plasma con- 
ditions and dynamics in accreting systems. 

We have presented our analysis of time-averaged grat- 
ing spectra of LMC X-4 in the high state and state tran- 
sition of its 30.3 day superorbital period. In both states, 
we detected a number of recombination emission lines; by 
comparing the intensities and ionization states of these 
lines at multiple epochs, we found evidence for two dis- 
tinct emission regions, which we identified as the X-ray 
photoionized stellar wind (Region I) and the outer accre- 
tion disk or stream (Region II) . We argued that the vari- 
ability of these regions, which was similar to variability 
observed in Her X-1 by JG02, was a natural consequence 
of accretion disk precession. 

We also found evidence for relativistic Doppler split- 
ting of the iron line in the inner accretion disk, and possi- 
ble splitting in low-Z recombination lines. We suggested 
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that these lines constitute a third region, the inner ac- 
cretion disk warp. Although the Doppler shifts should 
be much harder to detect in the high state, the variation 
of the equivalent width of th e iron line complex with ^E* 
(as measured bv lNaik fc Paull (|2003l) ) is consistent with 
our conclusions. 

While it is plausible that we simply observed two dis- 
tinct phenomena in the different observations, the pre- 
cession of the accretion disk provides a simple and coher- 
ent explanation for all the phenomena that we observed. 
As the disk precesses from edge-on to face-on, the con- 
tinuum flux rises; emission lines in the accretion flow 
brighten accordingly. At this point, the neutron star is 
visible and pulse-phased spectroscopy reveals the rapid 
motion of the accretion column, dragged along by the 
pulsar's magnetic field. As the disk turns its edge to our 
line of sight, the neutron star is obscured, but the ra- 
dial velocity of the corotating disk is maximized, and we 
observe Doppler splitting in the disk emission lines. 

In this work, we have performed a new high-resolution 
study of X-ray emission line variability with superorbital 
phase in LMC X-4. Accordingly, we are able to report 
a new spectral confirmation of accretion disk precession 
as the origin of the superorbital period. The variable 
Doppler shifts from the inner accretion disk and the rel- 
ative evolution of other emission lines in LMC X-4 are 
fully consistent with and predicted by the precession of 
the accretion disk. Future high spectral resolution ob- 
servations at the peak of the high state and during the 
low state will go a long way towards a systematic test of 
this model for accretion dynamics in X-ray binaries. 
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